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ABSTRACT 
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1 Introduction 



The sine-Gordon model is completely integrable (exactly solvable) both on the classical and 
on the quantum level We shall write the sine-Gordon equation in the following form: 

^^(^'^) ~ ^"(^'^) + ^sin/3M(a;,t) = . (1.1) 
Here m is a mass, (5 is the coupling constant. For later convenience we also introduce 

In the classical case u{x,t) is an function of two variables, x and t are space and time coor- 
dinates. In the quantum case u{x,t) is a local quantum field of the sine-Gordon model. The 
Hamiltonian reads 

n = ldx (^1/ + ^{d^ur + ^(1 - COS Pu)^ . (1.2) 

Momentum and topological charge are given by 

P = — J dx p dxU , Q = J dx dxU . (1.3) 

Here p{x,t) = dtu{x,t) and u{x,t) satisfy Poisson brackets {p{x),u{y)} = 6{x — y). Equation 
( |1.1| ) has a Lax representation and a classical r-matrix lip-0. After quantization, the fields u 
and p satisfy canonical commutation relations [M(x),p(y)] = i 5{x — y). The physical ground 
state \VL) of the quantum system can be obtained by filling the Dirac sea of negative energy 
pseudoparticles 0. 

Let us now consider the quantum operator 

exp(a(5(x)) =exp|^(M(x) -m(0))1 , Q(x) = ^ [ dzd^u{z) , (1.4) 

[ ZTT J ZTT JO 

where Q{x) measures the topological charge on the interval [0,x]. In this paper we show how 
to represent the correlation function 

{n\exp{aQ{x))\n) (1.5) 

as the determinant of an integral operator (in fact we shall see below, that the coefficient a in 
(|1.5|) needs to be renormalized) . Note that via differentiation with respect to a we can obtain 
correlation functions of local quantum fields from ( |1.5| ). We shall consider the quantum version 
of (|1.2|) in the region f < 7 < ^ (many of our intermediate results hold in larger regions of 
coupling constant). Note that 7 — is the quasiclassical region of the sine-Gordon model 
and at 7 = 7r/2 the spectrum of the Hamiltonian is equivalent to free fermions. To deal with 
the ultraviolet divergences of the continuum model we shall employ a suitably chosen lattice 
regular ization. 

The determinant representation then permits to describe the correlation functions in terms 
of a system of integrable integro-differential equations. These equations can be solved by means 
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of a Riemann-Hilbert problem which in turn enables one to obtain elementary formulas for the 
asymptotics of the correlation functions. This program has first been applied to the nonlinear 
Schrodinger equation in |§] and is described in detail in the book (see also jlOl ). 

There has been previous work on determining correlation functions in the sine-Gordon 



model. Form factors were determined by Smirnov in [0, |12[. At the free fermionic point 



7 = 7r/2 a determinant representation of the correlation function ( |1.5| ) has been constructed 
using the coordinate Bethe Ansatz in [jl3|. A description of a different correlator at the free 
fermionic point through a Fredholm determinant (derived from a form factor sum) which in 



turn satisfies an integrable differential (sinh-Gordon) equation has been obtained in [14|. In this 



paper we start the investigation of correlation functions in the sine-Gordon model for general 



7, in particular away from the free fermionic point in the framework of its solution [T^ by 
means of the Quantum Inverse Scattering Method. 

The plan of this paper is as follows: in section ^ we review the integrable lattice regulariza- 
tion of the sine-Gordon model introduced in The Algebraic Bethe Ansatz is formulated 
and the construction of the ground state |T^ is discussed. In section ^ we derive the deter- 
minant representation of the correlator ( p..5|) for the range of coupling constants stated above. 
As this part of the analysis is very similar to the analogous problem for the spin-| Heisenberg 
XXZ model (which was treated in full detail in |2^) we omit many details and only give an 
account of the main steps without providing prrofs (which can be found in [^). In sections ^ 
and 1^ we embed the determinant representation into a system of integrable integro-differential 
equations and in section ^ the related Riemann-Hilbert problem is formulated and the leading 
asymptotic behaviour of the correlation function is extracted. 



2 Lattice Sine— Gordon 
2.1 ^-Operator 

We shall consider a lattice version of the sine-Gordon model which is also completely integrable. 
It will have exactly the same r-matrix (both in the classical and quantum case) as the continuous 



model. The elementary /^-operator of the LSG model is |15, ^ 



g-i/3p„/8 -i/3p„/8 smh(X - i(3un/2) 

Cin\\) =1 I (2.1) 

-imA sinh(A + i(3un/2) e'^vn/ii p^^^iPPn/s 

Here A ist the lattice constant and p„, Un are the dynamical variables on site n of the lattice. In 
the quantum model they obey canonical commutation relations [M„,Pm] = i^nm- Furthermore, 
we have introduced 

1 ^2 



Pn = (I + 2S COS Pu^Y , 5=(^-mAj . (2.2) 

The symmetries of the /^-operator of the LSG model are expressed by the identities (the asterisk 
means Hermitian conjugation of the quantum operators) 

ay C*{n\\) ay = C{n\X) , a' C{n\\) a' = C{n\\ + n^) . (2.3) 

Its quantum determinant [|l^, |l^ is 

deiqC{n\\) = 1 + 2^cosh2A . (2.4) 
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The >C-operator ( |2.1| ) satisfies the Yang-Baxter equation 

R{X, /i) iC{n\X) ® Cinlfi)) = (^(nl/i) ® C{n\X)) R{X, /i) . (2.5) 
R{X, /i) in Eq. ( |2.5| ) is the standard sine-Gordon i?-matrix given by the following expression: 

\ 



i?(A,/i) 



//(/i,A) 

^(/i,A) 1 

1 gifi,X) 

V f{fi,X)J 



(2.6) 



Here 



/(/i,A) 



sinh(/i — A — 27) 



^(/^, A) 



sin7 



(2.7) 



sinh(/i — A) sinh(yU — A) 

In different sites of the lattice the matrix elements of C commute. As usual in the Quantum 
Inverse Scattering method (QISM) we define the monodromy matrix by taking products of the 
/^-operators in matrix space: 



r(A) = C{L\X) C{L - 1|A) ■ ■ ■ i:(l|A) 



(2.8) 



where L is the number of sites in the lattice which we take to be even. By construction this 
operator also satisfies a Yang Baxter equation 



i?(A, /i) (T(A) ® T{fx)) = (T(/x) ® T(A)) i?(A, /i) . 



(2.9) 



It might be interesting to point out that the entries fo the ^-operator ( p.l|) form a repre- 
sentation of a quantum group: The operators (we suppress the site index n) 



s- 



i sin 7 m A 

-2 
i sin 7 m A 

g-i/3«/2 ^1 



gj/3p/8 p gj/3p/8 



g-i/3p/8 p g-i/3p/8 



3j/3«/2 



satisfy the commutation relations of the quadratic (Sklyanin) algebra 





= ' ( 


S J s 


q — q ^ ^ 





^°)2 - (S 



1\2 



with q = exp(z7). For q being a root of unity this algebra has finite dimensional cyclic repre- 
sentations: for rational values of the parameter j/ir = Q/P the quantum operators entering 
the /^-operator can be written as 2P x 2P matrices with elements 



X 



Sate 



m{a-l)/P 



TX 



5„ 



a, h 



,2P ,a + 2P = a 



The definition of the C operator alone does not determine a definite lattice model: In 
addition the Hamiltonian of the lattice sine-Gordon model needs to be specified. For this 
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choice there exist several different possibilities (see |T^, |T6|, |T^). All of them are completely 
integrable and can in fact be diagonalized simultaneously. Furthermore all of them have the 
same continuum limit ( |1.2D . They differ from one another by higher orders in the lattice spacing 
A. While all of them can be considered equivalently as a lattice regularization of the continuum 
model we shall show below, how a unique lattice Hamiltonian can be chosen by requiring that 
it has the "same" ground state wave function as the continuum model. This choice of the 
Hamiltonian will bring the dynamics of the lattice model as close as possible to that of the 
continuum model. 



2.2 Algebraic Bethe Ansatz for the lattice sine— Gordon model 

We shall consider the monodromy matrix (|2.8| ) 

As a direct consequence of the Yang-Baxter equation ( p.5|) for T(A) the trace of the monodromy 
matrix, the so-called transfer matrix 

r(A) = trace T(A) = A(A) + /^(A) (2.11) 

commutes for different values of the spectral parameter A, i.e. [t(A), T{fi)] = 0. Hence, it is the 
generator of commuting integrals for the system which are diagonalized by the algebraic Bethe 
Ansatz. Starting point is the "pseudo vacuum" (or reference state): To construct this simple 
eigenstate of r(A) we combine the ^-operators in pairs: 

C{n\X) = C{2n\X) C{2n - 1\X) ^ ( fj^^ ) . (2.12) 



Choosing 



{u\0)n= U ~2ScOS-{u2n + U2n-l)> 5 ( U2„ - U2n-1 " - + ^ j (2.13) 



(for rational 7/7? = Q/P the (5-function can be replaced by a Kronecker 5-symbol and |0)„ will 
become normalizable) we find from (|2.12|) 



7n(A) |0)„ = 

«n(A) |0)„ = {l + 2^cosh(2A-27)}|0)„ (2.14) 
5n(A) |0)„ = {l + 25cosh(2A + z7)}|0)„ 

Now we can follow the standard steps of the algebraic Bethe Ansatz. As a consequence of 
( |2.14| ) the "global pseudo vacuum" 

L/2 

|0) = n I0)n (2.15) 

n=l 

is an eigenstate of the operators ^(A) and D{X) (and hence the transfer matrix ( |2.11| )) with 
eigenvalues a(A) and d{X), respectively: 

a(A) = {1 + 2^cosh(2A -27)}^ , rf(A) = {1 + 2^cosh (2A + ^7)}* . (2.16) 
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More eigenf unctions of the transfer matrix are found by acting with the operator -B(A) on the 
pseudo vacuum 

N 

UBiXm (2.17) 

provided that the {Xj} satisfy the Bethe Ansatz equations 

/ l + 2gcosh(2A, -27) \ * _ _ fr sinh(A,- - A^ + i-f) 

\1 + 2S cosh(2Aj- + Z7) / /J^ sinh(Aj - A^ - ^7) ^ ' ^ 

The corresponding eigenvalue of the transfer matrix (|2.11|) is 

N N 

A(A|A,) = a(A) n /(A, A,) + d{\) [] /(A„ A) (2.19) 

where /(A,/i) has been defined in ( p.7| ). 

The number of the Bethe Ansatz roots Aj can be identified with the topological charge 
( |1.3| ). The correct lattice version in the quantum case is 

^ L/2 ^ _ 

^ = 7^ E («2n - M2„-l) + L^^—^ . (2.20) 
n=l ^7 



The difference in the coefficient compared to (1^) is related to the fractional charge of the 



excitations. In it was shown that the fractional charge appears due to the repulsion 
beyond the cutoff in the process of ultraviolet renormalization. ( |2.20|) is the number operator 
for particles 

TV N 

Ql[B{\,)\0) = Nl[Bi\,)\0) . 
One can prove that (here is the Pauli matrix in the matrix space) 

Now we can discuss our choice of the lattice Hamiltonian for the lattice sine-Gordon model. 
As mentioned above we want to construct a lattice version resembling the dynamics of the 
continuum model as closely as possible. According to the standard quantization of the sine- 
Gordon model the ground state of the continuum model contains no bound states (strings). 

For possible lattice models we shall concentrate on the two integrable models introduced in 
Refs. [|1^, |T6[ and ||18|. The latter has been constructed by Tarasov, Takhtajan and Faddeev 
(TTF) such that it contains interactions of nearest neighbours on the lattice only. The ground 
state for this Hamiltonian was found in Ref. EO]: in addition to a Dirac sea of elementary 



particles it contains bound states. In the continuum limit the density of the bound states 
vanishes, thus reproducing the known results for the continuum model. Apart from the Hamil- 
tonian the QISM yields higher integrals of motion. These describe interactions over larger 
distances. Adding these interaction terms to the TTF Hamiltonian with coefficients vanishing 
in the continuum limit A — > produces different lattice Hamiltonians with the same continuum 
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limit while preserving integrability. This is the origin of the freedom in choice of the lattice 
hamiltonian. 

Another Hamiltonian for the lattice sine-Gordon model has been introduced in [l^, |T^: 



The corresponding ground state for this Hamiltonian has been constructed by Bogoliubov |T7 
He was able to prove that in the interval 7r/3 < 7 < 27r/3 the ground state is built from 
elementary particles only — ^just as in the continuum model. Furthermore, he found that the set 
of observable excitations coincides with the continuum model. Hence, unlike the situation in 
the TTF model no phase transition is met in performing the continuum limit. For the reasons 
stated above we choose this Hamiltonian for our studies of correlation functions. 

It is given in terms of trace identities. Expressing the zeroes d{K±) = and a{i'±) = of 
(gJ§ as 

2 9 

g2«± _ _^±ie-n ^2u± ^ _5±ien ^^ere b = , (2.21) 

(A± = I {i7T ± In 6) are the zeroes of the quantum determinant (|2.4|) of C) the Hamiltonian of 
the lattice sine-Gordon model considered here is given by 

326sin7] [dX a{X) 



dX d{X) 




+e--(|rln^l -e-(^ In -)-:() } (2.22) 



This is the model studied in [|T5|, |T^. From ( p.l9| ) one finds that ( p.lTp are eigenf unctions of 
this Hamiltonian with energy eigenvalues given by 

N / N \ N 

nLscUBi^m = EMA.) n^(^.)io) (2-23) 

j=l \k=l I j=l 

with the single particle energies 

^^^^ _ m^A I e^^ _ e-^T 

326i I sinh(K+ — A) sinh(fi;_(_ — Xj — 27) sinh(/t_ — A) sinh(K_ — A — 27) 

+ . (2.24) 

sinh(z/_|_ — A) sinh(z/+ — A + ^7) sinh(z/_ — Aj) sinh(z/_ — A + 27) J 

In the continuum limit A — (which is reached by letting 6 — here) one immediately 
reproduces the result 

^(A)|a^o = -m^Asin7 cosh2A 

for the single particle dispersion of the continuum model. 

To find the solution of (|2.18|) corresponding to the ground state of the model it is necessary 
to classify the possible configurations of Xj in the complex plane according to the so called 
string hypothesis The details of this are not important in the present context. It was 

found by Bogoliubov [O that the ground state of (|2.22|) is obtained by filling all permitted 
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states of pseudoparticles with rapidities Xj on the hne ImA = tt/2. Taking the logarithm of 
Eq. ( p.l8| ) in such a state one obtains 

L TT, ^ /sinhfAo — Ajfc + ^7) \ 

-p(A,+,-)^2.Q,-,i:i-'( ,i„4,; J ■ (2.25) 

Here the Qj are distinct integers characterizing the state uniquely and 

TT, , /l-25cosh(2A-i7)\ 
^^^ + ^2) = -^^Hl-25cosh(2A + .7)j ' 

In the thermodynamic limit the density p{Xj) = j^dQj/dXj is then given in terms of the integral 
equation 

-p'{X + Z-) = 27rp(A) + dfiK{X - /x)p(/x) (2.26) 

where 

K(X) = -sin 27 ^ -2 sin 27 

sinh(A + ^7) sinh(A — ^7) cosh 2A — cos 27 
This integral equation can be solved by Fourier transform resulting in (A = — In b) 

1 r ,1. .-ikx sinh |fc(7r - 7) 1 



^ Att J -00 sinh ^k'y cosh |/c(7r — 7) 2 ^ ^ 

Similarly one can compute the excitation energies. This is useful to find the correct mass 
renormalization formula: To perform the continuum limit of the sine-Gordon model one should 
let A — i> and simultaneously m ^ 00 as 



m 



const. A-^/^ . (2.29) 



3 Algebraic formulation of correlation functions 



For the evaluation of the correlation function ( |1.5|) we shall make extensive use of the similarity 



(in the framework of the QISM) of the LSG model with the spin-i XXZ Heisenberg chain which 
is derived from a monodromy matrix satisfying a Yang-Baxter equation with the same i?-matrix 
as the present model. The correlation functions corresponding to ( p.. 51 ) in the XXZ model 



have have recently been studied in |22, 



First we note, that the symmetry of the £ operator 



implies for the ground state 



configuration consisting of rapidities {Xj = Xj + z7r/2} with real Xj 



{bCx, 



C(A,) 



In order to express the correlation function ( |1 . 5| ) in the algebraic framework outlined above we 
first need to define the lattice analog of the operator Q{x) in ( [1.4|) . The correct expression is 
found to be 

^l(^) = ^ H (^2fc - M2fc-l) +^ I ^. 1 (3.1) 



k=l 



27 
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which counts the number of particles in the interval [1, n] {n even). In the continuum limit this 
expression becomes 

Qiin)^^{u{x)-u{0)) + ^(^^^^ , x = nA (3.2) 
Hence the lattice analog of the correlation function ( [L.5D can be written as 

/oi f nf ^Mn^ (0| nf=i CjX,) exp{aQ^{n)) ULi B{h)\0) 

(iZ expfaOWn)) iZ) = r; (3.3) 

^ ' ^ ^ ' (oinf=iC(A,)nf=i5(A.)|o) 

where Xj are solutions of the Bethe Ansatz equations ( p.l8|) for the ground state configuration. 

Let us first study the norm appearing in the denominator of this expression. To evaluate this 
expression one should commute the C{Xj) to the right of the product where they annihilate the 
pseudovacuum |0). Since the commutation relations between the elements of the monodromy 
matrix (|2.10|) are completely determined by the i?-matrix we can use the result of [0, (see 
also [ p5| , p6| ) for the norm of Bethe Ansatz states (after identifying 7 with 2{7c — rj) in paper 
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iV TV I I f ^ _ _ I 

(0| n C(A,) n B{Xkm = i-sm^f n /(A.' n «(A,)rf(A,) detA^ (3.4) 

i=i J=l [j^k J [j=l J 

where the N x N matrix A/" is given by 

J^jk = S,, ln^ + J2 K{X, - A„)| - K{X, - A,) . 

{ dXj d[Xj) n=i J 

The functions A'(A) and ci(A), d{X) have been introduced in the previous section. In the 
thermodynamic limit this expression can be further simplified: We rewrite M = 1 - J where 

_ K{X, - Afe) _ 

-^jr'fc "jk Q J Ujk ^jk^j 

^^ = '^^^1k\ + T.k{x,-k). 

OAj a[Aj) n=l 

Comparing the last expression with Eqs. ( p. 251) and (|2.26| ) for the ground state density of 
particles one obtains 6j = —2%Lp{Xj). Performing the thermodynamic limit on the matrix X 
one finds that it turns into a Fredholm integral operator X = 1 + acting as 

X * /U = /(A) + — / diiK{X - p)/(/x) (3.5) 

Here A'(A) is kernel given in (|2.27| ). 
Putting everything together we find 

{Q\X{C(X,)XlB{\m = (27rLsin7)^ IIM'A' 

X < n «(^.)^(^.)P(A.)| det (1 + • (3.6) 
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We now turn to the numerator of (|3l3|): to reduce the evaluation of the expectation value of 
exp{aQi{n)) in a Bethe state ( p.l7| ) to the computation of scalar products we divide the lattice 
of length L into two sub-chains of length n and L — n and associate a monodromy matrix with 
each of them, namely 

r(A) = T(2, A)r(i, A) , r(z, a) = ) , z = i, 2 . (3.7) 

In terms of /^-operators they are given by 

T(2, A) = C{L, A) C{L - 1, A) ... C{n + 1, A) 
r(l, A) = C{n, A) £(72 - 1, A) ... £(1, A) . 

By construction these monodromy matrices satisfy the same Yang-Baxter equation ( p.9|) as 
T(A). Similarly, the global reference state ( p.l5| ) can be decomposed into a direct product of 
pseudo vacua for the subchains |0)2 ® |0)i (remember that we have chosen n to be even) which 
are eigenstates of Ai{X) and D.i{X) 

A(A)|0)i = a,(A)|0), , A,(A)|0)i = d,(A)|0)i, (3.8) 



where ai(A) and di{X) are given by ( 2.16 ) with L replaced by n and L—n for i = 1,2, respectively. 
The creation and annihilation operators -Bj(A) and Cj(A) act according to 

a(A)|0), = 0, (0|5,(A) = 0. (3.9) 

In this decomposed quantum space the numerator of (p.3|) can be rewritten as (see e.g. 

E i(oi nci(A?,) n5i(Afjio)i 2(oi n c^ixna) n ^^(Af.jio)^ 

Ic Ib lie IIb 



Jb,Ic ) [iIbJIc 

n /(Af.,Af,j|( n (3.10) 

Jb,IIb ) [icllc ) 

where the sum is over all partitions 

{Afj U fXfjJ = {A}, {Xfj n {Af,J = , {Afj U {A?,J = {A}, fxfj n {Af,J = 

of the set {A} with cardjA/^} = cardjA/^} = rii, card{A//(^} = cardjA/Zg} = N — rii. Due to 
( |3.9| ) we only need to consider partitions such that the sizes of Ib and Ic (and IIb and lie) are 
the same. We next turn to an investigation of the scalar products occurring in ( |3.1CI|) . Owing to 
(|3.8|) and (|3.9| ) and the fact that the monodromy matrices T(i, A) fulfill the same Yang-Baxter 
equation ( |2.9| ) as T(A) it is sufficient to consider scalar products on the entire lattice 

N N 

SN={0\l[C{Xf)l[B{Xm ■ 

j=i k=i 
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Here we do not assume that the sets of spectral parameters {A^} and {A*"} are the same, and 
we also do not impose the Bethe equations (p.l8|) . From ( |2.9|) and the action on the reference 



state y4(A)|0) = a(A)|0), D{X)\0) = d{X)\0) it follows that scalar products can be represented 
as 



A,Dj=l k=l J 1^ // 

where the sum is over all partitions of {A'"}U{A^} into two sets {A^} and {A^}. The coefficients 
Kn are functions of the Xj and are completely determined by the intertwining relation ( p^.9| ). 
The i?-matrix ( |2.6|) is however identical to the one for the spin— | Heisenberg XXZ model 
(after appropriate identifications of the coupling constants). This implies that the coefficients 
K]\f for the sine-Gordon model and the XXZ chain are identical, so that we can take over the 
result for the XXZ case (see e.g. [^]). The main point is that the K^^s can be represented 
as determinants. This is done in two steps: first the so-called highest coefEcients, which are 
obtained for the partition {A"^} = {A*-^}, {A^} = {A^}, are represented as determinants 

*"({A°} {A«}) = ('n9(\-^f)9(^f.\°)jnMAf,Af)det(M|), 



where from 



sinh(A - fi-h) , . -sin'^-f 
/i(A,/i) = — , t(A,/i) 



—i sin 7 ' ' sinh(A — — i'y) sinh(A — /i) 

In the second step arbitrary coefficients Kn are then expressed in terms of highest coefficients 
as follows 

^4\xn Ia^!) = fn n/(Af,Ar)Un n /(a^a^) 

Vl^ J 1^ ;/ \jeACkeDC J \l&ABmeDB 

/{A^-} {A-^}W^ ({X^^} {A--}, .3^3. 



Using ( p.l2|) and ( p.l3| ) in (|3.11|) we obtain the following expression for general scalar products 



in the lattice sine-Gordon model 

= n^?(Af,A^)^7(Af,Af)Esgn(Pc;)sgn(PB)nMAf ,Ar)nMAr,Ar) 

j>k j,k l,m 

X n MA^ , AD n MAf , A^ )det(M^^)det(M^g) , (3.14) 

l,k j,m 

where Pq is the permutation {A^*", . . . , X^'^\ Af", . . . , A^^„} of {Af , . . . , A^}, Pb is the permu- 
tation {Af ^, . . . , A^^, A^^, . . . , A^^„} of {Af , . . . , A^}, sgn(P) is the sign of the permutation 
P, and 

(m^^) = tiXf^X^diXnaiXf). (3.15) 
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Following the steps first carried out in ||2^ it is now possible to represent Sn as a single 



determinant. The discussion for sine-Gordon is identical to the only for the XXZ chain [22 



so that we only present a brief discussion of the necessary steps and give the final result. We 
first note that the sum on the r.h.s. in ( 3.14 ) looks very similar to a Laplace decomposition 
of the determinant of the sum of two matrices {Si)jk = t{\'j'\ X^)a{Xf)d{\^) and {S2)jk = 
t(Af , A^)(i(A^)a(Af ) (see e.g. p. 221). However this does not reproduce the /i(A, /i)-factors. 
This leads to the introduction of a dual quantum Geld acting in a bosonic Fock space with 
vacua |0) and (0|^ according to 



y.(A) = MA) + g(A), [y,(A),y.(/x)]=0, (0|g(A) = = p(A)|0) 
[p(A),g(/x)] = -ln{h{X,^i)h{^i,X)) , [p{X),p{fi)] = 0=[q{X),q{fi)] . 



(3.16) 



We emphasize that the field ip commutes for different values of spectral parameters. Using the 
dual field it is now possible to recast ( |3.14| ) as a single determinant of the sum of two matrices 



N 



s 



N 



n 9{Xf, A^)^(Af , Af ) n a{Xf)diXf) J] hiXf , X 

j>k j=l j,k 



det5|0) 



S. 



•jk 



t(Ay,An+t(A^Ay)^exp(^(A 



r(Af 

''iiiMA^,Af)MAf,A^)' 



(3.17) 



where r(A) 



The price we pay for representing Sn as a single determinant is the 



occurrence of the expectation value in the dual space. 

Using (|3.17|) in (|3.1CI|) and then applying the dual field trick several times it is possible to 
represent (p. 10) as a single determinant of the sum of four matrices. This analysis is completely 



analogous to the XXZ case treated in [22| so that we only state the result: 



TV 



(oi nc^(A,)exp(«Qi(n)) n B{Xkm = n /(a., a 



k=l 




N 



na(A,)rf(A,)^(0|detg|0) 



Gjk = t{Xj, Afc) + t{Xk, Xj) ^^^-^^ exp (v?2(Afc) - v^2(Aj)) 



riiXk 

+ exp (a + if^iXk) - V53(A 
d 



^(Afc, Xj) + t{Xj, Xk 



ri(Aj 
riiXk 



exp 



[(fi{Xj) - </?i(Afe)) 



/ 



—% 5jfcsin7 



dX, 



N 



ln(r(A,)) + Eln 



V 



n=l 



h{Xj, Xn) 



(3.18) 



h{Xn, Xj] 

\ ^ and the commuting dual fields are defined 
^a{X) = Pa(A) + ga(A) , (0|g,(A) = = pa(A)|0) ,(0|0) = 1 , a = 1...4 , 



where r'i(A) = ai{X) / di{X) 
according to 



/ l+25cosh(2A-t7) \ ' 
Vl+25cosh(2A+j7), 



^We use the same notation as in fE2| 
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/I 1 ON^ 

10 1 

111 

VI 1 1/ 



ln(/i(A,/i)) + 



/I 1\ 
110 
10 11 

Vo 1 1 ly 



(3.19) 



where a, 6 = 1 ... 4. Here all terms not proportional to 6jk in Qj^ are understood in the sense of 
I'Hospital for the diagonal elements. In the thermodynamic limit further simplifications take 
place. Following the analysis for the norms above we express Q as the product of two matrices 
J and W 

g = -{sm-f)WJ , Jjk = ^jkOk , Wjfc = 6,k - ^V(A,-, A^) , (3.20) 



where 6j = ~2tx Lp{\j) and 



(sin7)V(A,/i) = t(A,/i) + t(/i,A) 



n(A) 



+exp (a + v24(/i) - V^3(A)) 



exp (<^2(/i) -^2(A)) 
r,(A) 



t(/i, A) +t(A,/i)— -exp ((^i(A) - V5i(/i)) 
n(/i) 



.(3.21) 



In the thermodynamic limit W turns into an integral operator W = 1 + acting as 

dfiV{\^i)f{^) , (3.22) 

where the integral kernel is obtained from (|3.21|) as (the arguments of the dual fields are shifted 
by in/ 2 which does not alter the defining commutation relations (|3.19| )) 



;i + ^v^)*/|a = /(a) + ^ 



r(A,/i) 



— sm7 



1 



+ 



e2^(A)e2(/i) 



sinh(A — /i) I sinh(A — /i — ^7) sinh(A — /i + ^7) 

1 er^(/i)ei(A) 



+ exp(a + v24(/i) - V^3(A)) 



+ 



with 



e2(A) 



1 -25cosh(2A + i7) 



3</'2(A) 



sinh(A — + ^7) sinh(A — /i — ^7) 
1 -2^cosh(2A-27) 



(3.23) 



,l-2Scosh(2A-z7), 
Putting everything together we thus find 

((]|exp(aQi(n))|(]) ; 



ei(A) 



l-2^cosh(2A + 27) 



(6|det(l + ^V)|0) 
det (1 + IK 



(3.24) 



where 1 + -^V and 1 + i^K are integral operators acting according to (|3.22|) and (^]5|) with 
kernels defined in (|2.27| ) and ( p.23| ). 
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4 Continuum linit 



As mentioned in the introduction the purpose of the present work is to determine correlators for 
the SG Quantum Field Theory, and the lattice model studied above is used merely as a regulator 
for the UV divergences. We are therefore interested in the continuum limit of the determinant 
representation (|3.24|) . As mentioned above the SG Quantum Field Theory is recovered from 



the lattice regularization by taking the lattice spacing to zero A — > and simultaneously the 



bare mass m to infinity keeping mA^ fixed |17]. In order to take the continuum limit we 



now employ the following regularisation for the integral operators in ( p. 241) : we restrict the 



integration for the integral operator 1 + to the interval [—A, A], and then take A — in 
such a way that S'cosh(2A) ^ 1 VA G [—A, A] (recall ( p.2|) for the relation of S and A). Using 
this regularisation the ej(A)'s simplify to 

e2(A) = exp(-ipsinh(2A) + <^2(A)) , ei(A) = exp(ipsinh(2A) + y?i(A)) , (4.1) 

where ^ 

p = ^La^ sin(7)nA . (4.2) 



Here we have used ( p.29| ) and nA = x should be identified with the continuum distance. The 
constant c is given in terms of the physical soliton mass. 

This regularisation allows to embed the determinant ( p.24| ) into a system of integrable 
integro-differential equations which we shall need later to determine the subleading terms in 
the asymptotic expansion of the correlation functions: with ( |4.1| ) the kernel ( p.23| ) can be 
brought into standard form 0]: We perform a change of variables z = exp(2A), and replace 
the factors (sinh(A — ± ^7))^^ in ( |3.23| ) by an integration over an exponential. Then the 
transpose of the kernel (|3.23|) reads (up to a similarity transform which leaves the determinant 
unchanged) 

—V^{z,,Z2) = / dsJ2Ej{z2\s)e,{z^\s) , (4.3) 

27r zi- Z2 Jo j^i 

where 

^^i(^l^) = ^^exp(v?4(2;))|2,^,s) , Ei{z\s) = y=exp{-ip3{z)){2, z, s\ 

v27r \/Z7i 



e2{z\s) = -^\l,z,s) , E2{z\s) = -^{l,z,s\ 
v Zir V Zn 



(^3{z\s) = -^exp{-ipk{z) + ^2{z))\2,z,s) , E3{z\s) = ^exp{ipk{z) - ip2{z)){2, z, s\ 

v27r y Zn 



eiizl-s) = —=exp{-ipk{z) - ipi{z) + ip4{z))\l,z,s) 
V Ztt 



Ea{z\s) = -^exp{ipk{z) + (pi{z) - (p3{z)){l,z,s\. (4.4) 

V ZTT 

Here we use the notation k{z) = \{z — z~^), w = exp(z7), k = exp(^), and 



\l,z,s) = ^ 2z sin(7) exp{izws) = {2,z,s\ , \2,z,s) = ^J2z sin(7) exp(— z— s) = {l,z,s\ 

(4.5) 
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are normalized in such a way that (1|1) = ds (1, z, s|l, z, s) = 1, and similarly (2|2) = 1. 
The inverse of the integral operator 1 + ^V'^ is defined by 



;i - ^) * (1 + i-r) = 1 = (1 + ) * (1 - i?) ^ 



271 



In terms of the functions fj{z\s), Fj{z\s) 



[l-R)* e, 



f,{z\s), E,*{l-R) 
the kernel of R can be written in a form similar to (|4.3|) 

Rizi,Z2) = / /i(^lk)-^i(^2|s 

^1 ~ ^2 .'-1 "'O 



FMs) 



as can be seen by acting with (1 + ir^^) (E3 



(4.6) 



(4.7) 



(4.8) 



5 Integro- Differential equations 

Let us now derive integro-differential equations determining the functions fj{z\s) and Fj{z\s). 
The analog of these equations in the case of impenetrable bosons proved very useful for the 
anlysis of the corresponding Riemann-Hilbert problem and we expect the equations below to 
play a similar role for the problem at hand. To this end we consider derivatives with respect 
to p and the integration boundary A. For the A-derivatives we find 



dAfj{z\s)+J2 dtU,i{z\s,t)fi{z\t) = 
1=1 •'^ 

dAFj{z\s)-Y. dtFi{z\t)Uij{z\t,s) = 
1=1 •'^ 



where 



Ujkiz\s,t) 



32A 



-2A 



z — z — e~ 

The p-derivatives of the functions fj{z\s) obey the integro-differential equations 

4 r 



dpfj{z\s) 



-ik{z)fj{z\s) + 
\ 1=1 



r{0) , 1 



{Sj,3 + 5i,4) 



Z,S I 

4 



k=-i 1=1 ^fc=3 



(5.1) 



/,(e-^^|s) F,{e~'^\t) . (5.2) 



(5.3) 
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where Ijk{s,t) = 6jkS{s — t) and where the integral operators B^"^^ and C'-"-* are defined as 

/ \ /•exp(2A) 

Ck\s,t) = / -f){z\s)E,{z\t) . (5.4) 

v/cxp^— 2A) Z 

We note the following relations between the integral operators Bj^^ and C]^"* 



(5,z5(s-t). (5.5) 



These identities can be easily proved by using ( [4.7| ). From now on we will replace B^^f} in all 



expressions by C]?. The IDE for F. j{z\s) are found to be 



dpFj{z\s) 



'tk{z)F,{z\s) + lj2Fi* 



1=1 



^(0) , 1^(1) 



-^^^:^:Fl*[I+^c^%k*B 

k=3 1=1 

The "potentials" B^"^ and C^'"'^ obey the equations 



2,S I 

4 



(5.6) 



( ) -'- (0) ( ) 



m=3 



Z, j^-^ 



mi 



C%-'\s.t)-C%-''\s,t)+^Y.Cf*C\l 



y{n+l) 



(0) . 



(5.7) 



m=3 



4 4 



s,t 2 



+ ■^{^,3 + ^,4) 



5 



m=3 i=l 

jk 



Im 



* B 



(1) 

mk 



s^t)-B%-''\s^t)-^Y.Bf,UC, 



{n) , r^{0) 
Ik 



1=1 



B 



'jk I 



.,t)-i?]r^H.,t) 



The derivatives with respect to A are given by 

dACf^is,t) = 2e'^f,{e'^\s)Fk{e'^\t) + 2e-'''fj{e-'^\s)Fk{e~'^\t) 

dAC^l\s,t) = 2/,(e2^|.)(^F,(e2^|t)+^Fz*Q^)[^^J+A^-A 



aACf(.,t) = 2e-^V,(e^l3) F,(e^^|t)+zF,*C« 



9Ai?«(.,t) = 2F,(e2^|t) /,(e2^|.)-^i?«*/, 



+ A ^ -A . 



(5. 



A ^ -A , 
(5.9) 
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Eqs. ( |5.1| ), ( |5l3| ) and ( ^.6| ) form a Lax pair. Their consistency is implied by the following 
relation for the cross-derivatives 



dpdAfj{z\s) = dAdpfj{z\s) . 

In order to simplify the computations we first introduce some notation. We rewrite 

i ^ 1 ^ 

dpfj{z\s) = --zfj{z\s){6j;s + 5i,4) + Y.^Ji* fi ^ji * fi 

^ l=\ ^ Z=l 

where 

In the same notation (|5.6| ) can be written as 



(5.10) 

as 

(5.11) 



dpFj{z\s) = -zFj{z\s){6j,3 + 6j,i) - ^ * aij 
^ 1=1 

Similarly we introduce the notation 

Ti r.k .,MMs^.M-A\s,t) 

Uji{Z\b, I) — ^ 



2A ' 



(5.12) 



(5.13) 



(5.14) 



z — e^'^ z — e 

where Ajk{A\s,t) = 2ie^^ fj{e^^\s) Ffc(e^^|t). In what follows we will denote by Ajk{A) the 
integral operator in the s-variable with kernel A{A\s,t). After some calculations we arrive at 
the following equations 



iz 



dAdpfj{z\s) = — (5j,3 + 5j,4) 



m=l 



4 ^(^)jm * fm 

E 



m=l 



z — e 



2A 



+ A ^ -A 



V 



J. /^A,„(A)*/„ 
1=1 



\m=l 



z — e 



2A 



A ^ -A 



+ 



m=l 



,4 / 4 

^ «=1 \m=l 



^im(A) * /„ 



z — e 



2A 



A -> -A 



(5.15) 



dpdAfj{z\s) 



E 

r?i=l 



apA(A)j>„ * /„ 



z — e 



2A 



^ m=l 



^(A)j,m * /„ 



z — e 



2A 



v 



+ A ^ -A 



+ A ^ -A 



(5m,3 + <5m,4) 



E E f ^^r% + A - -a) * («z™ * + -&z™ * 

i=\ m=i \^ e y \ Z 



(5.16) 
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In order to equate ( p.l5| ) and ( p.l6| ) we first rewrite both equations in the form Ojm * fm, 
where O are complicated integral operators, and then "truncate" the /^'s from the resulting 
expressions, which amounts to supposing that they form an independent set of functions in the 
space the integral operators act in. In the next step we then compare the resulting expressions 
(which are both meromorphic functions of z) at the singular points z = oo, 0, e^^"^. If they 
(their residues) are equal at these points the expressions coincide for all values of z. For z ^ oo 
we get the condition 



dAajm{s,t) = -{Sm,3 + ^m,4 " ^j,3 " Sj^4){Aj^{A\s,t) + Ajm{-A\s,t)) 



'2 A 



1=1 



A^-A 



At 2; = we obtain 

At z = we get 

apAj„(A|s,t) 

+ E (K-« + e^'^^^id * ^im(A) - Aji{A) * [ai^ + e-^\ 



(5.17) 



(5.18) 



-je^^{5m,Z + 5mA - Sj,3 - 4) Aj>„(A | S, t) 



lin 



(5.19) 



s,t 



whereas the condition from z = e~^^ is obtained by taking A —A in (^.19 ). It is straightfor- 
ward to show that these equations hold by inserting the expressions for a, b and A and using 
the identities for the p- and A-derivatives of Cj^'' and 5^^^ written above. 

Finally, to relate the functional determinant in ( |3.24| ) to the quantities introduced above 
we turn to the logarithmic derivatives of det(l + ^V'^)- 

The derivative with respect to p is given by 



a^ln det(l + — 



tT{il-R)*—d,V' 
Zn 



(5.20) 



Using 



we find that 
1 



27r 



dpV {zi,Z2) 



k{zi)-k{z2) 



Zl - Z2 



ds ej{zi\s)Ej{z2\s) 

4 



1 1 r°° 

-(1 + )/ ds J2e,{z,\s)E^{z2\s). (5.21) 

Z Z1Z2 Jo 



This implies that 
(1 - ^) * ^d,V 



1 

2 Jo 



ds ^ 



exp(2A) 



j=3"'exp(-2A) 

+ 2Jq II/j(^iM)^i(^2|s). 



dz [S{zi — z) — R{zi, z)] 



ej{z\s)Ej{z2\s) 



ZZ2 



(5.22) 



i=3 
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Using the representation ( [4.8| ) of R{zi,Z2) we rewrite the r.h.s. as 
r.h.s = - / ds y dz[d{Zi — z) — R{Zi, z)\ — 

2 Jo ,„ Jcxv(-2A) 



n poo poo 

- ds 
2 Jo Jo 



Jcxp(-2A) Z1Z2 
k=3 1=1 -'cxp{-2A) 



Z2 



\ roc 4 

n ds J2fj{Zi\s)Ej{z2\s) . 

Z JO ^-^3 



Using this with ( ^.4|) in ( 5.201 ) we finally arrive at 



a,in(det(i+^^-))4|:/;.. 



.(2) 



(1) M^) 



Ik 



1=1 



The logarithmic derivative of the determinant with respect to A is 



9Aln(det(l + — y 
2n 



(5.23) 



(5.24) 



(5.25) 



After some manipulations similar to the case of the Bose gas (see [0) this can be rewritten as 

-2A 4 ^ 



-4A 



R{e 



2 sinh(2A) 



-2A „-2A\ 



E / ds fj{e 
.7=1 



Y: / dt /Ke2^|t)FKe-2^|t) 



2A 4 



U=l 



(e-'^^ls) 



,4A 



2 sinh(2A) 



^ fOO ^ /-oo 

E/ rf. /,(e^^|.)F,(e-^^|s) dt fiie-'^mie'^^lt) 

7=1 -^0 Li=i 



(5.26) 



This embeds the determinant into the system of integrable integro-differential equations derived 
above. 



6 The Riemann-Hilbert Problem 

In this section we show that the results of the previous section can be reformulated in terms of an 
infinite-dimensional Riemann-Hilbert problem (RHP) for an integral operator valued function 
Y[z). This connection will enable us to determine the asymptotic behaviour of the correlation 
function. We introduce the conjugation matrix G{z) of this Riemann-Hilbert problem as 



[G{z\s,t)]^^ = 6ij6{s -t) + 27re,{z\s)Ej{z\t) . 

It's elements can be expressed in terms of the projectors ( |4.5| ), e.g. 

[G{z\s,t)]^-^ = 5{s -t) - K'^exp{(fi{z) - Lp3{z))\2,z,s){2,z,t\ 
[Giz\s,t)]^^ = Kexp{(p4{z))\2,z,s){l,z,t\ 



(6.1) 
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Consider now an integral-operator valued function Y{z) with kernel Yjk{z\s, t), j, k = 1, . . . ,4, 
s, t G [0, oo) acting on a vector / of functions of z and s according to 

/■oo ^ 

[Y{z)*f{z)]^= dtY^Y^,{z\s,t)Mz\t) . (6.2) 

k=i 

Y{z) is solution to the following RHP 

. r(^)=/ + Er=ii^forz^oo. 

• Y{z) is analytic throughout the complex plane with the exception of the contour C, which 
is the interval [exp(— 2A), exp(2A)] on the real axis (see Fig. [l|). 

• Y'{z) = Y^{z)G{z) on C where Y^{z) are the boundary values as indicated in Fig. ^ 
and ^(2;) is the conjugation matrix ( |6.1|) . 

This RHP can be rewitten as the system of singular integral equations 



as 



The solution of ( |6.3D can be expressed in terms of the functions E and / defined in Section ^ 

pM2A) ^_^ J^{z'\s)E^{z'\i 

/exp(-2A) Z' — Z — iO 

which follows from the identity 

f,{z\s) = / dt Yjk{z\s,t) ej{z\t) . (6.5) 



The potentials B'-^^ and (|^) can be related to the solution Y{z) of the RHP through 
asymptotic expansions around and 00. We find 

Yjk{z) — > I + tC^^'> +tzC^^'> +tz^C^^^ + 0{z^) foTz^O, (6.6) 
Yjk{z) — >I- + 0{z~^) for z^oo. (6.7) 



From (|6.6| ) and ( ^.5| ) we find 

[/_,5(i)]*C'(^) = _,(y-i(0)A 

dz 



2=0 



Y{z)) . (6.8) 



Together with ( p.24| ) this expresses the correlation function ( ^.23| ) in terms of the solution Y{z) 
of our RHP. 
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6.1 Analysis of the RHP 



While the operator- valued RHP defined above determines the correlation functions completely, 
its solution appears to be a daunting task in general. In what follows we concentrate on the 
leading term in the asymptotical decomposition of the solution of the RHP in the region of 
coupling constant f < 7 < The reason for this restriction is the following: the upper bound 
on 7 stems from the construction of the ground state of our lattice regularization. The lower 



bound ensures that the parameter p defined in ( |4.2| ) will go to infinity in the continuum limit, 
which essentially simplifies the analysis of the RHP: it permits us to study the asymptotical 
decomposition of the solution of the RHP with respect to p (recall that p contains the continuum 
distance as well). Due to the fact that this parameter will be not only large but diverge the 
number of terms in the asymptotical decomposition will be very small - in fact we expect 
only three contributions (see also below). As we shall show in our analysis of the leading 
contribution, the special form of the conjugation matrix in addition to our interest in partial 
traces of Y only allows to reduce the RHP to a tractable scalar one (still containing the auxiliary 
dual fields, of course). The analysis of the subleading terms is technically much more involved 
and is currently under investigation. We will report on this work elsewhere. 

Let us now turn to the calculation of the leading term. First, we note that the conjugation 
matrix can be decomposed into the product of an upper and lower triagonal matrix as follows 



Here 



^ roo 
c=l 



/I ai{z\s,t) a2{z\s^t) eyiY>{ipk{z)) a^{z\s,t) eyi-p{ipk{z))\ 



(6.9) 



Ti{z\s,t) 






VO 



1 ai{z\s,t) eyiY){ipk{z)) a^{z\s,t) eyi]i{ipk{z)) 
1 aQ{z\s,t) 

1 



(6.10) 



with matrix elements 

ai(z|s, t) 
a2{z\s,t) 
a:i{z\s,t) 
a4^{z\s, t) 
ar^{z\s,t) 

aQ{z\s, t) 



Kexp((y94(z)) 



1 + ^2 exp(v94(z) - (f'iiz)) 



|2,z,s)(l,2;,t| 



-- G^^{-^2{z) + ^^{z))\2, z, s){2, z,t\ , 



K 



exp{(pi{z) - (p3{z) + (Piiz)) 



2,z, s){l,z,t\ 



1 + ^2 exp{ip4{z) - (psiz)) 
— exp{-ip2{z) + (p3iz) - ip4{z))\l,z,s){2,z,t\ 



K 



/texp(v9i(z) -y^3(^)) ^ 
1 + hi'^ exp{ip4^{z) - ip3{z)) 
Kexp{^2{z) + ^i{z) - ^3{z)) ^ ^ 
1 + k'^ exp{(p4{z) - ip3{z)) 



(6.11) 
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Similarly, we find 
/ 

T2{z\s,t) = 



ci{z\s,t) \ 

(3iiz\s,t) C2{z\s,t) 

P2{z\s,t) exp{—ipk{z)) P4{z\s,t) exp{—ipk{z)) C3{z\s,t) 

\P3{z\s,t) exp{—ipk{z)) P5{z\s,t) exp{—ipk{z)) PQ{z\s,t) C4{z\s,t) J 



(6.12) 



The matrix elements of T2 are given by 



1 



and finally 



Pi{z\s,t) 


= 


P2{z\s,t) 


= 




= 


p4{z\s,t) 


= 


P5{z\s,t) 




(3(i{z\s,t) 




Ci{z 


s,t) 




s,t) 




s,t) 


Ci{z 


s,t) 


go throng 


a a ' 



-exp{-(pi{z))\l,z, s){2,z,t\ 



Kexp{ip2{z) - V23(^)) 



-\2,z,s){2,z,t\ 



1 + k"^ exp{(p4{z) - (fsiz)) 
-K'^exp{-Lpi{z) - Lpsi^z) + Lpi{z))\l,z,s){2,z,t\ 
exp{ip2{z)) 



■\2,z,s){l,z,t\ 



1 + k'^ exp{(p4{z) - (p3{z)) 
Kexp{-ipi{z) + ip4{z))\l,z,s){l,z,t\ , 
-Kexp{-ipi{z) -'f2{z) + ip4{z))\l, z, s){2,z,t\ 



(6.13) 



1 - 

1 + 
1 - 



|2)(2| 



k"^ exp{Lp4^{z) - 03(^)) 
1 + k"^ exp{Lp4{z) - ^siz)) 

\exp{-Mz)+V3iz))\l){l\ 
1 



1 + ^2 exp{ip4{z) - (p3{z)) 
1 + expi^^iz) - ^3{z))\l){l\ 



|2)(2| 



(6.14) 



We define an integral-operator valued function Y{z) in the following way: 

• Y{z) = Y{z) outside the "bubble" defined in Fig 2. In particular Y{z) = Y{z) for 
z 0, cxD, which will be important later. 

• Y{z) = Y{z)Ti[z) in the region enclosed by the real axis and the contour Fi. Note that 
in this region lmk{z) > Vz. 

• Y{z) = Y{z)[T2{z)]~^ in the region enclosed by the real axis and the contour F2. Note 
that in this region lmk{z) < Vz. 

It can be easily seen that the function Y{z) defined in the above way has the following prop- 
erties: Y{z) is analytic in the whole complex plane with the exception of the contours Fi and 
F2. On the contours Tj Y satisfies the conjugation equations 



(Ynz) 
(Ynz) 



YHzmz) 

Y^{z)T2{z) 



zeT,, 
zeT2 . 



(6.15) 
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Since we are only interested in the asymptotic behaviour of the determinant for p ^ 1 we can 
use the fact that in this hmit Ti(2) become blockdiagonal in the vicinity of the contour ri(2) 
from which we find that ~ 

yi^)-rf^ . (6.16) 



$2(^) 

Here ^j{z) are solutions to 2 x 2 operator- valued RHPs 

^{z)=^{z)*G,{z) , J = 1,2 (6.17) 
with the same conjugation contour C as the original RHP and conjugation matrices 
^ \2){2\ \ ( ^B^|2)(l| \ 

^H^J ^ 1 - ^ l^_£iPkM£))|l)(2| (l + 2SP(£3(£te(£)))|l)(l| j ' ^^■^'^) 



G2{z) 



1-|2)(2| 

i-|i)(i| 

Kexp{ipi{z) + ip2{z) - ip3{z))\2){l\ 

-Kexp{-(^l{z) - Lp2iz) + Lp4{z))\l){2\ (1 + k"^ exp{(^i{z) - V93(^)))|l)(l| 



Using the fact that Gj{z) form representations of GL{2\C) we can now calculate the deter- 
minants of Gj{z) as is shown in the appendix 

det(Gi(2;)) = exp(-a + Lp3{z) - (fiiz)) , det(G'2(2)) = exp(a - Lps{z) + (p^iz)) . (6.19) 

The scalar RHPs for the determinants 

det($-(^)) = det($+(z))det(G,(z)) , j = 1,2 

is now easily integrated to give 

1 .exp{2A) _a + ipsiz^) - ipi{Zi) . 



det($i(2;)) = exp(-— / dzi 

2m Jexn(-2A) 



exp(-2A) Zi — Z 

det($2(2:)) = exp(-— / dzi ) . (6.20) 

2m Jexp(-2A) Zi — Z 

7 Leading Term in the Asymptotics of the Correlator 

Let us now relate the solution of the scalar RHPs to the logarithmic derivative of det(l + ^^"^)- 
The contribution due to C-"^ in ( p.24| ) can be obtained from ( |6.7| ) and ( |6.2CI| ) as 



lJ^r°°, ^fnu . iz 

I /.exp{2A) 
An Jexp(-2A) 



i:J^^dsCS^{s,s) = jim|ln(det($2(^)); 

dz[a - ips^z) + (p4{z)] . (7.1) 
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Similarly, the second contribution in ( p.24| ) is with ( |6.7| ) 

2dz 



ln[det($2(^))] 

2 = 



AtT Jcxp{-2A) 2;' 



Combining these to the leading asymptotical bahaviour of (9pln(det(l + ^v)) and using the 
fact that they are p-independent we obtain 

1 /. 1 rrT^ A ,Ci;psinh(2A), ,p /■exp(2A) 1 X/ / X /XXX 

det l + — = Aexp -^^ ^ exp^/ dz {I + -){^^{z) - ^-.{z))) 

/TT TT 47r Jcxp{-2A) Z'^ 



= ^exp(^^^^^^^^^)exp(^ r dAcosh(2A)(v94(A) -(^3(A))) , (7.3) 
vr vr J-A 

where A is a p-independent constant and where in the last step we have changed back to 
the original A-variables. Decomposing the combination of dual fields into "momenta" and 
"coordinates" and using the commutation relations ( |3.19| ) we find 

^4(A) - ^3(A) = P(A) + g(A) , [Qif,), P{\)] = . (7.4) 

This enables us to trivially evaluate the expectation value with respect to the dual fields in 
this approximation: the dual fields are found not to contribute at all leading to the following 
result for the leading asymptotical behaviour of the correlator 

{n\exp{aQi{n))\n) ~ Iexp(— sinh(2A)) , (7.5) 

TT 

where A is a constant independent on p. 

We will now argue that the approximation ( [7.5[ ) is too crude due to the fact that we have 
neglected the influence of the dual fields in the subleading factors in the solution of the RHP. 
We expect the final answer for the solution of the RHP to be of the form 

det(l + i-\/^) = C({^,}) exp(C({^,})ln(p)) 

X exp |^ Q^Ps^^M2A) ^ j\xcosh{2X){ip^{X) - ^3(A))^ ,(7.6) 



where we keep in mind that p — * 00 as the lattice spacing A — * 0. In ( |7.6|) C is p-independent 
and we have conjectured that the subleading term in the solution of the RHP is a power-law 
in p. Evaluating the expectation value of (|7.6|) in the dual bosonic Fock space the dual Gelds 
will contribute in the exponential term, i.e 

{exp{aQi{n))) ~ (o| det(l + ^t/^) |0) 

= C p^ exp {jnp) exp (^pln(p)) exp — ^ — ^ + cjjap^ . (7.7) 

Here ^ and m are functions of 7, the soliton mass etc. For this answer to be of the correct 
qualitative form, the following conditions have to be satisfied: 
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^ = 0, as the leading asymptotic behaviour should be exp(const p). 

The last factor in ( |7.7| ) has to be cancelled by a suitable regularization procedure for the 
result to make sense. In the continuum limit we have (|3.2|) which implies that 



(exp(agi(n))) ^ (exp (^j[u{x) - «(0)])) x exp (^^ . (7.8) 

We see that this expression contains a divergent factor depending both on a and on the 
distance x. We now adjust our "cutoff" A in such a way that the divergent factor in ( |7.7|) 
precisely reproduces the divergent factor in ([r.8|), i.e. 



faxTr — j\ / sinh(2A) \ 

exp [-^^ ) = exp (^[^^ + • 

If a; = this leads to the following relation between the "cutoff" A and the lattice spacing 
A 

exp(2A) = A-'^ , \ . (7.9) 
with a finite constant c. The procedure outlined above fixes the relation between A and 



A and thus between the divergent part of p and A. Note, however, that the result ( |7.9| ) 
for this relation is not consistent with the requirement 5'cosh(2A) -C 1, which we have 
used in order to simplify the kernel of V in section ^. Therefore the assumption u = 
has to be wrong and we do need a A-dependence of u instead which corrects (|7.9|). 



8 Summary and Conclusion 

In this paper we have applied the method of to correlation fucntions in the sine-Gordon 
model. In order to deal with the UV divergences we used an integrable lattice regulariza- 
tion of the sine-Gordon model to derive a determinant representation for quantum correla- 
tion functions. We then took the continuum limit and obtained a determinant representation 
for the sine-GordonQFT. Furthermore we embedded the determinant in a system of inte- 
grable integro-differential equations which we showed to be associated with an operator-valued 
Riemann-Hilbert problem. The quantum correlation function was expressed in terms of the 
solution of this RHP. We then presented a general approach to obtain the leading asymptotical 
behaviour of the solution of the RHP, which in turn yields the leading term in the asymptotics 
of the quantum correlation function. We showed that the subleading terms in the asymptotical 
decomposition are essential for obtaining explicit expressions for the asymptotics of the correla- 
tion function due to the presence of the dual quantum fields. For the case at hand there appear 
to be only two subleading terms in the asymptotical decomposition which is very encouraging! 
The analysis of the subleading terms is a difficult mathematical problem by itself and we will 
report on it in a separate publication. 
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A Gl{2\C) Representation by Integral Operators 

In this appendix we show how to essentially simplify the analysis of the operator-valued RHP 



through the use of GL{2\C) representation theory. We closely follow the discussion of |^ 
Let us consider an integral-operator valued 2x2 matrix with kernel 



0{s,t) 



s,t E [0, oo). 



(A.l) 



02l{s,t) 022{s,t) 

Multiplication of integral-operator valued matrices O and V is defined in the usual way as 

2 



[OV],^is,t) = Y. dr 0,kis,r)Vkjir,t) , t,j = 1,2. 



(A.2) 



The left (right) action of the integral operators Oij on fuctions defined on the interval [0, oo) 
is given by 



dt Oijis,t) fit) , g*0,j 



t JO 



ds f{s) Oij{s,t) . 



(A.3) 



Let us now construct a special class of such operators O which form a representation of Gl{2, C): 
we start with two pairs of functions {a{s), l3{s)) and {A{s), B{s)) on [0, oo) which we represent 
in Dirac notation as a{s) = |1), /3{s) = |2), A{s) = (1| and B{s) = (2|. These functions are 
chosen in such a way that 

POO 

(1|1) = / dsA{s)a{s) = 1 = (2|2) = 



dsB{s)l3{s) 



In this notation we may write left multiplication by Oik as 



dik\i) 



dt0^k{s,t)a{t). 



(A.4) 



(A.5) 



Observe now that one may define a representation A of Gl{2, C) in terms of integral operators 
via 



Me G/(2,C) I — >AiM) 







/-|2)(2| 





( Mn 




Mi2|l)(2| \ 












1 M21 


2)(1| 


M22|2)(2| ; 



(A.6) 



Here Mu, M12, M21 and M22 are complex numbers and / is the identity operator in the space 
of integral oprators on [0, 00). Multiphcation by the integral operators |1)(1|, |1)(2|, |2)(1| and 
|2)(2| is given by e.g 

\l){2\f{s)=^fdsB{s)f{s))\l). (A.7) 
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Therefore act like projectors on the "states" \i) and 

In particular identities like [/ — |1)(1|] |1)(1| = are seen to hold. Indeed for any M, e 
Gl{2, C) the representation A has the following properties 

(PI) A{MN) = A{M)A{N) ; A{I) = I ; A{M~^) = ^^^(M) 



(P2) 



/ 



Tr 



A{M) 




trM = Mil + M 
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(P3) 



Det^(M) = det M = M11M22 - M12M21 . 



(Ai 



Properties (PI) and (P2) can be established by direct computation using the rules given above. 
Property (P3) shows that the determinant of the integral operator A is simply equal to the 
determinant of the 2x2 matrix M, which is quite remarkable. It is established by expressing 
the determinant as a trace via InDet^ = Trln^, then using (PI) in the expansion of the 
logarithm, using (P2) to express the operator trace in terms of the matrix trace, and finally 
expressing the sum over traces back as determinant of the matrix M. 

It can be easily checked that the representation ( |6.18D of the conjugation matrices Gj(z) is 
precisely of the above form (here z plays the role of a parameter), which in turn allows us to 
evaluate the determinants of the conjugation matrices. 
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